//.do file for Effects of teaching experience and culture on choral directors’ descriptions of choral tone//
// Frizzell, Emily Y. and Windsor, Leah. PLOS-ONE 2021.//
// For inquiries: elyoung1@memphis.edu//

// use Frizzell-Windsor-PLOS-ONE.csv //

********************************************************************************
//Base model regression results for Figures 1-5, and 7-8
********************************************************************************
********************************************************************************
ologit healthy wc auxverb negate adj compare posemo negemo discrep as sex2 agegrp experience2 ethnicityrace, ro cl(id)
eststo
ologit appropriate wc auxverb  negate adj compare posemo negemo discrep as sex2 agegrp experience2 ethnicityrace, ro cl(id)
eststo
estout, cells("b(star label(Coef.)) se(label(Std. err.))") stats(r2 N, labels(R-squared "N. of cases"))  label legend varlabels(_cons Constant)
estimates clear
********************************************************************************

********************************************************************************
//Topic model regressions for Figure 6
//AS1 Contemporary Western Choral
//AS2 Asian subcontinent (Indian/Bollywood)
//AS3 3:1 (Spiritual, American/Contemporary, Jazz/Folk) //different version of familiar song
//AS4 Western/European Choral
//AS5 South African Musical
//AS6 Bulgarian Traditional
//AS7 American Folk, Traditional
//AS8 Black Gospel
********************************************************************************
********************************************************************************
ologit healthy wc auxverb  negate adj compare posemo negemo discrep as sex2 agegrp experience2 ethnicityrace, ro cl(id)
margins, at(as==1) at(as==2) at(as==3) at(as==4) at(as==5) at(as==6) at(as==7) at(as==8)
marginsplot, yline(0)

ologit appropriate wc auxverb  negate adj compare posemo negemo discrep as sex2 agegrp experience2 ethnicityrace, ro cl(id)
margins, at(as==1) at(as==2) at(as==3) at(as==4) at(as==5) at(as==6) at(as==7) at(as==8)
marginsplot, yline(0)
********************************************************************************
